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^ ; Abstract 

This paper is concerned with theoretical analysis of a heat and moisture transfer model 
arising from textile industries, which is described by a degenerate and strongly coupled 
parabolic system. We prove the global (in time) existence of weak solution by constructing 
an approximate solution with some standard smoothing. The proof is based on the physcial 
nature of gas convection, in which the heat (energy) flux in convection is determined by the 
p I ' mass flux in convection. 



< 



< 



Key words: Heat and moisture transfer, Porous media, Global weak solution. 

1 Introduction 



Mathemaitcal modeling for heat and moisture transport with phase change in porous textile 
^ . materials was studied by many authors, e.g., see O [6l [9l [12]. A typical application of these 

I models is a clothing assembly, consisting of a thick porous fibrous batting sandwiched by two 

thin fabrics. The outside cover of the assembly is exposed to a cold environment with fixed 
temperature and relative humidity while the inside cover is exposed to a mixture of air and 
, vapor at higher temperature and relative humidity. In general, the physical process can be 

' viewed as a multiphase and single (or multi) component flow. In this process, the water vapor 

■ moves through the clothing assembly by convection which is induced by the pressure gradient. 

The heat is transferred by conduction in all phases (liquid, fiber and gas) and convection in gas. 
Phase changes occur in the form of evaporation/condensation and/or sublimation. Based on 
^ ' the conservation of mass and energy and the neglect of the water influence, the model can be 

described by 



^(ea) + ^(nea) = -r,e, (1.1) 



I {eC.,MC.T + (1 - e)C.sT) + ^ {eC^.MuC.T) = ^ 



Here C„ is the vapor concentration (mol/m'^), T is the temperature (K), e the porosity of the 
fiber, M the molecular weight of water and A the latent heat of evaporation/condensation in 
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the wet zone while in frozen zone, it represents the latent heat of sublimation. C^g and C^s are 
the heat capacities of the gas and mixture solid, respectively. 

The evaporation/condensation (molar) rate of phase change per unit volume is defined by 
the Hertz-Knudsen equation [10] 



Rf V 2ttRM V VT 

where R is the universal gas constant, Rf the radius of fibre and E is the nondimensional 
phase change coefficient. The vapor pressure is given by P = RCyT because of the ideal gases' 
assumption. The saturation pressure Psat is determined from experimental measurements, see 
Figure ??. 

The vapor velocity (volumetric discharge) is given by the Darcy's law 

kkrg dP . , 

u = ^— (1.4) 

[Ig OX 

where k is the permeability, k^-g and pg are the relative permeability and the viscosity of the 
vapor, respectively. 

Numerical methods and simulations for the heat and moisture transport in porous textile 
materials have been studied by many authors with various applications [H [U [HI [T71 [19]. How- 
ever, no theoretical analysis has been explored for the above system of nonlinear equations. A 
simple heat and moisture model was studied in [T8], where the model was described by a pure 
diffusion process (without convection and condensation) with a non-symmetric parabolic part. 
There are several related porous media flow problems from other physical applications. A pop- 
ular one is a compressible (or incompressible) flow in porous media with applications in oil and 
underground water industries, which is described by an elliptic pressure equation coupled with 
a parabolic concentration equation for incompressible case and a system of parabolic equations 
for compressible case. The existence of weak solution for incompressible and compressible flows 
has been studied in [3l[7| and [1], respectively. However, in most of these works, the temperature 
is ignored and the phase change (condensation/evaporation) does not occur due to the nature of 
these applications, while both temperature and phase change play important roles in the textile 
model. 

For the textile model, the water content in the batting area usually is relative small and one 
often assumes that all these physical parameters involved in the system (|l.ip - (ll.2p are positive 
constants. With nondimensionalization, the system (jl.ip - (|1.2p reduces to 

Pt - {{pO)xP)x = -F, , . 

{pe)t + aOt - {{pOhpOh - (^^x)x = AF, , ^ 

for X G (0,1), t > 0, where (•)^ = ioi p = x,t, p = p{x,t) and 9 = 9{x,t) represent the 
density of vapor and the temperature, respectively. 



il/2 



Ps{0) 



and PsiO) ~ Psa,t{0) /O^/'^ . a and A are given positive constants and k = k\ ^ K2(? is the 
heat conductivity coefficient with (i = 1,2) being positive constants. We consider a class of 
commonly used Robin type boundary conditions [Sj [HJ [H [20] defined by 

{pe),pU=i = a\p' - p{l,t)), (A/>U=o = aO(/>(0,t)-/>0), (1.6) 
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and 

ne^[,=i = p\e^ -e{i,t)), Ke.U=o = /?°(e(o,t) -^O), (1.7) 

and the initial condition is 

p(x,0) =po(x), e{x,{))=eo{x), xe(0,l) (1.8) 

where a^,a^ represent the mass transfer coefficients, p^,p^ are the density of the gas in the 
inner background and outer background, respectively, P^,(3^ the heat transfer coefficients, and 
9^,6^ the inner and outer background temperatures. Physically, all the parameters above are 
positive constants and po{x) > p,0o{x) > 6, where p and are positive constants. Based on 
the experimental data in Figure ??, we assume that ps is a smooth, increasing and nonnegative 
function defined on M"^ which satisfies 

lim — - — = 0, lim — — = oo (1.9) 

for some > 0. For physical reasons, we set Ps{Q) = for < 0. 

The objective of this paper is to establish the global existence of weak solution to the initial- 
boundary value problem (jl.5p - (|1.8p under the general physical hypotheses (|1.9p for the saturation 
pressure function F. The difficulty lies on the strong nonlinearity and the coupling of equations. 
To the best of our knowledge, there are no theoretical results for the underlying model. More 
important is its significant applications in textile industries. Also analysis presented in this 
paper may provide a fundamental tool for theoretical analysis of existing numerical methods. 
Our proof is based on the equivalence of mass and heat transfer in convection. 

2 The main result 

Before we present our main result, we introduce some notations. Let T be a given positive 
number in the following sections. We define 

O = (0,l), /=(0,r], Qj = J7x(0,t], Qt = 5^x/, 
Vi{Qt) = L\l-H\n)), V2{Qt) = {/ G L\Qt) | ||/||y,(Q^) < +oo}, 

II/IIv2(Qt) = ess sup ||/(-,t)||L2(C) + II/x||l2(Qt)' 

te[o,T] 

wI'\Qt) = {/ G L\Qt) I ftJxJxx G L\Qt)]. 

Let 2?(i7 X [0, r)) be the subspace of C°°(M^) consisting of functions which have compact support 
in M X [-oo,r), restricted to H x [0,r). 

Now we give the definition of weak solution to the system (jl.Sp - lll.Sp and then, state our 
main result. 

Definition 2.1 (Weak solution) We say that the measurable function pair (p, 9) defined on 
X [0, T) is a global weak solution to lll.5]) - [T78\) if {p,9) £ {Vi{Qt))'^ and the density p and the 
temperature 9 are nonnegative functions satisfying 

T f-T 

a'^{p{0,t) - p°)(f>{0,t)dt+ / a''{p{l,t) - p^)<j3{l,t)dt 
Jo 

+ / {-p(l)t + {p9)xp(t>x + T^)dxdt = / po(l)odx (2.1) 
Jo Jn Jq 
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and 



[ [a°(p(0, t) - p°)0(O, t) + f3\e{0, t) - ^°)]V(o, t)dt 



Jo 



Jo 




+ / / - (P^ + + {pO)xpOipx + kO^iP^ - XTtl; dxdt 



' {poOo + (y0o)'>Podx 



Jn '- 



(2.2) 



/or any test functions "0 G x [0, T)). 

Theorem 2.1 // the initial value {po,Oq) satisfies po E L^'^'^{Q,) (7 > 0), 6*0 S and 
Po ^ 0, ^0 ^ ^ for some positive constant 6, then there exists a global weak solution {p,0), in 
the sense of Definition \2. i\ to the initial-boundary value problem l^l.5\) - [T73\) such that 



In the following sections, we denote by Cpj^pj^... ^p^. a generic positive constant, which depends 
solely upon pi,P2,- " ^Pk, the physical parameters Ki,K2,cr and A and the parameters involved 
in initial and boundary conditions. In addition, we denote hy C{pi,p2, • " iPk) a generic positive 
function, dependent upon the physical parameters ki,K2,o" and A and the parameters involved 
in boundary conditions, which is bounded when Pi,P2) ■ " " ^Pk are bounded. 

3 Construction of approximate solutions 

Throughout this section, we let e be a fixed positive number which satisfies 



and {) < u < £. To prove the existence of global weak solutions to the system ()1.5p - (jl.8p . we 
introduce a regularized approximate system as follows: 



plnpG L~(0,r;Li(17)), p£L\Qt). Px & L\Qt); 
6,6-' eL^iQr), {I + p)6, e l\Qt). 



0<e< m.m{^,p\e°,6\l}, 



Pt - {{e + {p9)y)p^)^ - {p{pe9x)e) 



PX%V9) + x%Ps{6)), 



{p9 + a9)t - {k'9^)^ - {{e + {p6),))pj% - {p{pe6x)e9)x 
= Xpx%V9) - Xx'iPsi9)) + iX + 9) ix'{Ps{6)) - , 



in Qt 



(3.1) 



where is a cut-off function defined by 




and 



At^ = Kl + K2\P£ 
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and the subscriptions e, v define the smoothing operators in general by = Ext(/) * with 
H = v,e. Here r/^ is the standard mollifier and Ext(-) is the extension operator which extends 
any measurable functions defined on to be zero on R^\fiT- 
The system (j3.ip can be rewritten as 

[pt- ((e + (p0).)p..)x. - {p{PeB^)e\ + PX'W~Q) = x'{Vsm. 

\ (p + a)et - - \{e + {pQ\)p^ + p{peQ.)e\ 0x - PxHV0)e + (A + e)ps{e) = Xpx'iVe) . 

(3.2) 

The corresponding initial and boundary conditions are given by 

(e + ipe)^)p^ + p{pe0cc)el^^ = a\p^ - p{l,t)), 

(e + ip0)u)p. + p{PeOcc)el^o = ""(/'(O'*) " Z^")' 

p{x,0) = poe{x) := {po)£{x) +e, 

(3.3) 

K'e,\,=i = pHe^-e{i,t)), 

'^"exU=o = /?°(^(0,t)-^°), 

e(x,0) =0o.(x) := {eo)e{x). 

We prove the existence of solutions to the system (|3.2p - ()3.3p by using the Leray-Schauder 
fixed point theorem. The following lemma (see [13], [H]) is useful in our proof. 

Lemma 3.1 (Aubin— Lions) Let Bi B2 ^ he reflective and separable Banach spaces. 
Then 

{u E LP (/; Si)| ut G L'?(/; B-,)} LP{I; B2), 1< p, g < 00; 

{u G L'?(/; B2) n Ll(/; Si)| G L1(/; B3)} LP{I; B2), l<p<q<oo. 

3.1 Existence of approximate solutions 

We define 

X = {ue L^{I; H^{n))\ -u > 0}, y = {n G wI'^{Qt)\ u > 0}. 

By Aubin-Lions lemma, Y X. Let e and be given positive constants and the parameter 
s G [0, 1]. For any given {p^,9^) G A^, we define p to be the solution of the following linear 
parabolic equation 

Pt - {{s + {p'e\)p,), - {p{pX)e). + spx'iVe^) = sx'ipsie'')), (3.4) 

with the initial and boundary conditions 

' {e + {pOe\)p, + p{pyje = aHsp'-p), at x = 1, 
{e + {p^e%)p^ + p{py,)s = a'^{p-sp^), at x = 0, (3.5) 
p{x,0) = spo£{x), for X G ri. 

Now with p in hand, we define 6 to be the solution of the linear parabolic equation 

{p + a)et - (K^e,), - [{e + {p°9\)p, + pipy)e] 

1_ (3.6) 

-spx'{Ve^)e + s{\ + 0)p,(0) = s\px'{^), 
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with the initial and boundary conditions 

' K^e^ = /3^{se^ -6), at x = l, 
n^e^ = p°{9 - s9°), at X = 0, (3.7) 
9{x,0) = s6oe{x), for x£Q. 
Let M denote the mapping from {p^,9^,s) to {p,9). Then we have the fohowing lemma. 

Lemma 3.2 The mapping M : X"^ x [0, 1] X"^ is well defined, continuous and compact. 

Proof. By the L^-theory of linear parabolic equations ^JJ, there exists a solution p G 
W2'^{Qt) for the system (I331)-(l33]) such that 

\\p\\w'^^^Q^)<C{e-\\\{p'9\\\c,^Q^^,\M^^^ 
By noting the fact 

\\pUm(n) < CeWpoWmn), Up'oWIc^q^^ < C.,T\\p'\\LHQr)P"\\LHQT)^ 

\\iP°e^x)e)x\\c^(Q^) < Ce,T\\p°&'^\\L\QT) ^ C'e.tIIp" ||l2(Qt) H"^!' IIl2(Qj,) , 

for the standard smoothing operator, we have 

Mw^-\Qr) ^ ^) (3-8) 

and therefore, 

Wph^iQr) < \\P\\,^2.^^Q^^ < C{e-\u-\\\p'^\\x,\\0°\\x,T). 

Let p'^ = max{p, 0}, = max{— p, 0}. Then p = p^ — p^ . By multiplying p^ on both sides 
of the equation (|3.4|) and integrating the resulting equation over Qt, we have 

/* [\e + {p^9\)\p-\^dxdT+ f [\sx'{V9)\p'\^ + sx'{Psmp')dxdT 
^ Jo Jo Jo Jo 

+ I (a°|p-(0,r)p + a°spV(0,r))dr+ [ {a^\p- {1,t)\'^ + ahp^ p- il,T))dT 
Jo Jo 

- / P' Px {p°eGx)sdxdT 

Jo Jo 
t fl 



lo Jo 

Notice that p^ > 0. Thus we have that 



1 II " '1 







n\p I'^dxdr. 



By Gronwall's inequality, we can see that p =0. Thus p = p^ > 0. This and (3.8) imply 
that p£Y X. 

Similarly, by the L^-theory of quasi- linear parabolic equations [11], there exists a solution 
9 G W^'^iQr) for the system (I3j])-(l321) and 



.2,1.^ ^<C{e-\u-\\\p'\\x,\\9'\\x,T). (3.9) 
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Let 9+ = max{(9,0}, 9' = max{-6',0}. Then 9 = 9+- 9'. Multiplying 9-/{p + a) on both 
sides of the equation (j3.6p and integrating the resulting equation over Qt, we can get 



Mr 



2 

t rl 



+ 
+ 
+ 




JO 

t L 



dx + 

sXpx'{V9^) 
-{0,T 



t rl 



t rl 



Jo P + ^ 



\9^ \^dxdT + 



JQ 



p + a 

It „ai 



-dxdr 



s{X + 9)ps{9) 
(p + a) 



9 dxdr 



p{l,T)+a 
ft rl 



[5\s9' + 9~{l,T))dT 



p(0,t) +0- 
t rl „ p- rt rl 

en- 




SPX% 



\9- 



-dxdr 



JQ 



K 



Jo 



-dxdr + 



JO 



p + a 

[{e + {p'9\)p, + p{pl9%]9- 



p + a 



-dxdr. 



Since Ps{9) = for 6* < 0, we observe that (A + 9)ps{9)9'' = a.e in 17-^. By Cauchy inequality 
and the estimations (3.8)-(3.9), we can estimate the terms in the right hand side of the above 
equality. Thus we obtain 

[' \9-\^dx <C{e-\u~\\\p'>\\xA\9'>\\x,T) f C\9~\^dxdT. 

JO JO JO 

Gronwall's inequality gives that 9^ = 0. Thus 9 = 9+ > Q. This and (3.9) imply that 9 G 
Y X. 

We conclude that the mapping M : X'^ x [0, 1] X"^ is a compact mapping. 
Now we prove the continuity of the mapping M. For any {p^,9^,s) G X"^ x [0,1], let 
{p,9) = M{p^,9^,s). Then 



(3.10) 
(3.11) 



pt - [{e + {p'9\)p, + pipy.hl + spx'iV90) = sx'{Ps{9')), 
{p + a)9t - {k'9,)^ - [{e + {p^9'^)u)Px + p{plel)e] 4 

-spx'{VT^)9 + s{\ + 9)ps{9) = rsXpx%^), 
with the initial and boundary conditions 

' (e + (pO^O),)p, + p(pO0O), = aMsy-;5), at x = 1, 

{e+{p^9\)p^ + p{p^A = a\p-sp^), at x = 0, 

p{x,0) = spos{x), for X £ Q. 



(3.12) 



and 



k'9^ - 



Denote p = p — p and 9 



l3^{s9^-9), at x = l, 
I3^{9 - s9°), at X = 0, 
9{x,0) = s9o!r{x), for x G . 

Then p satisfies the following equation, 

-x^(V^)] 



(3.13) 



where 



Pt-iF- F), + spx'W9^) + {s- s)px%V90) + sp[x% 
= is- s)x'{Ps{9')) + s[x'{Ps{9')) - xHPs{9'>))], 



F={e+ip"9").)p. + p{pyje, 



(3.14) 
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F = {e + {p'e'),)p^ + p{p^,e%. 

Multiplying the equation (j3.14p by p and integrating over Qt gives 



t)dx + 



t fi 



^0 



p. 



Iclxdr < C 



f f p^dxdr + (s - sf + 11/5° - pYx + - 0' 
Jo Jo 



0||2 
X 



With C = Cie-\i^-\ Wpoeh^n), ||/5°||x, ||^°||x,r). 

Thus Gronwall inequality implies that 



\\prx<c[{s-sr + \\p^-pTx + 

Similarly, we can derive the equation for 6 and get 



IX 



< C 



{s - sY + \\p'- 



a0||2 

\x 



+ 



f-eYx 



30 fl0||2 



Thus, the mapping M : X x [0, 1] X is continuous. The proof of Lemma 3.2 is complete. 



In addition, for s = we can see that M{p, 6*, 0) = for any (p, 9) G X"^ . Thus, by the Leray- 
Schauder fixed point theorem, there exists a fixed point for the mapping M(-, •, 1) : X"^ X"^ if 
all the functions (/?, 9) E satisfying 



{p,9)=M{p,9,s) 



(3.15) 



for some s G [0, 1] are uniformly bounded in X"^ . In fact, by the proof of Lemma 13.21 M 
maps {p,9,s) G X'^ x [0,1] into Y'^. Therefore, if {p,9) is a fixed point of M(-,-,l), then 
{p,9)eW^^\QT). 

Theorem 3.1 Under the assumptions of Theorem \2.1\ the system i f^.^j -i TOI) has a (strong) 
solution {p,9) G W2'^{Qt) which satisfies 

P>P^j, and 9t<9<9t for {x,t)eQT- (3.16) 

l|p||L°°(/;L4(n))) \\Px\\l'2{Qt), \\PPx\\l'2(Qt) < Cs,T, 

II^I|i,°°(Qt)' \\^x\\l'2{Qt), \\PeOx\\L'2{QT) < Ct ■ (3-17) 

where p and C^^t are positive constants which depend on e and T, independent of v; 9j', and 
9'r and Ct are positive constants, dependent upon T and independent of e and v. 

By the Leray-Schauder fixed point theorem, it suffices to prove the uniform boundedness of 
functions (p, ^) G X^ satisfying the equation (j3.15p and (|3.16p . 

3.2 Uniform estimates 

We assume that (p, 9) G and therefore, (p, 9) = M{p, 9, s) G Y'^, for s G [0, 1], i.e., {p, 9) is a 
(strong) solution of the following system, 

pt - ((e + {p9),)px)x - {p{pe9x)e)x + spx%V9) = sx%Ps{0)), (3.18) 

(p + a)9t - {k9x)x - [{e + {p9)u)px + p{peOx)e\ Ox 

- spx'iV9)9 + six + 9)psi9) = sXpx'iVd), (3.19) 
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with the initial and boundary conditions 

{e + {p9)u)px + p{Pedx)e = a^isp^ - p), at x = 1, 

{£ + {pe)^)px+p{peex)e = cP{p-Sp^), at X = 0, (3.20) 

p{x,Q) = spQe{x)-, for X G f], 

and 

K^^^ = /3i(s0i-0), at x = l, 

KfOx = P^{9 - 30°), at X = 0, (3.21) 

9{x,0) = s9o£{x), for x £ Cl, 

In this subsection, we derive some uniform estimates for solutions to the above initial- 
boundary value problems. 

Firstly we add the equation (jS.lSp multiplying by (A + 6) into (j3.19p and then, integrate the 
resulting equation over Qt. We arrive at 





x=l ,. 


f H2{x,t) 


dT < 


Jo 


x=0 Jo 



(It < {Xpoe + P0e0o6 + <T6'oe)(x)dx 
Jo 



where 



H2{x,t) = [epx + {p9)^px + p{pe9x)e]{^ + 0) + k'0x ■ 

With boundary conditions in (|3.2U|) - (|3.21|) . we have 

x=l 



H2{x,t) 



x=0 



= aHpil, t) - sp'){X + 6(1, r)) + a\p{0, r) - sff){X + 9(0, r)) 

+f3\9{l, t) - s9^) + (3\9{0, t) - s9°) 
> -a^sp^9{l, t) - a^sp^9{0, r) - Xsia^p^ + a^f) - s{p,W + p^9^) 



and therefore, 



[ {Xp + p9 + a9){x,t)dx <Ct + C [ 
Jo Jo 



')\\c(n)dT, 



(3.22) 



where 



Ct = {X + ||eo.||L-)||po.||Li + ^ll^oelli- + [A(aV' + aV°) + (Z?'^"' + /3°^°)] T. 



Similarly, subtracting the equation (j3.19p multiplied by 9^/1 from the equation (|3.18p multi- 
plied by 9^^^ /{I + 1) and integrating the resulting equation over Qt, we arrive at 



+ 







f H^{x,t) 
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x=Q Jo Jo 



+s{l + l) f [\x + 9)ps{9)9'dxdT= !\p^, + a){9oe)'^\x)da 
Jo Jo Jo 



(3.23) 



t rl 



JO 



19^+^ + X{1 + 1)9^] px%V9)dxdT + 



t rl 



x'{Ps{9))9^~^^dxdT. 



Jo 
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where 

x=0 

+ {l + l)PH0a,r) - se'Ml,T)f + + l)/3°(e(0,r) - s0"°)[0(O, r)]' 
= [aV(l,T) + (/ + - ahp^][eil,T)f+^ -{1 + l)phe^[9il,T)f 
+ [aV(0,r) + (/ + - aOsp°][0(O,r)]'+i - (/ + l)P%eyiO,T)]\ 

when / is large enough. Since 6^ < 9^/"^ + 6^+^ for any 61 > and /> 1, by (|3:22]) - (l3:23D . 

/ {p + a)e^+\x,t)dx + l{l + l) / / K^e'^^^pdxdr + s/ / / {\ + e)ps{e)eUxdT 
Jo Jo Jo Jo Jo 

<Ci^T + Cos [ [ l{l + e^+^'^)pedxdT 
Jo Jo 

t 



<CIt + CosI I \\e{;T)\t+2{ldT (3.24) 



lL°°(n) 

iJ u 

where 



1 



Jo 

and = C/^T + CI. RecaU the Gaghardo-Nirenberg inequahty 

<c'||/||i2(f,) + c'||/||^//(^)||/,||^/'(^), yf€H\n). 

i+i 

With f = 6 2 in the above inequahty, we obtain 



21+3 



and by Holder's inequality, 

* 21+3 21+3 

<ci f [\e'^)W.dxdT + j—^—— f [\\{e'-^u^dxdT 

Jo Jo + J^j<--o<--i Jo Jo 

/■* /"l (i + l)(2i+3) / + 1 ft 

<Cl / e^^+^ dxdT + / K^'-^^pdxdr. 
Jo Jo 4CoCi Jo Jo 

It follows that 

l\\0{-.r)tl'Ldr <^ f fe^^'I'dxdr + Cl f f o"^^ dxdr 
Jo Jo Jo Jo Jo 

+17^ f f ^iO'-'\e^A'dxdT. 

40o Jo Jo 
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By the assumption (T^ . we observe that CoC26''+i < ^^(0)6'^+^ + C for ah 61 > 0. Substituting 
the last inequahty into ()3.24p gives 



Ps{9)e'+'dxdT 



Jo Jo Jo ^ Jo Jo 

/■* (i+l)(2i+3) 

<C'iT + C3sl^ / e dxdr, (3.25) 

Jo Jo 

for / being large enough. Let Iq be a positive integer satisfying 
(/o + l)(2/o + 3) ,^.^2Vf2 

where is defined in ()1.9I) . By noting the fact 

(io + l)(2io+3) 1 

with li = 2{Iq + 2 + r])/r], we have 

Jo ^ Jo Jo 4 Jo Jo 

where C'l^ j> = C'l^ j> + CriClof^ for some constant Ct independent of /q- Furthermore, 

sup / e^''+\x,t)dx+ [ [ \{e'-^)^\^dxdt <C'i'^T 

0<t<TJo Jo Jo 

and by the Sobolev embedding inequality, 

rT 







Since Iq is a fixed positive integer dependent solely upon r], we obtain the estimate 

rT 



From (f3:22l) and (IHIM . we get 



0<t<T 

and 



/ ||0||Loc(f^)dxdt < Ct. (3.26) 

JO 

sup [ ip + pB)dx < Ct (3.27) 

<t<T Jo 



rl rt rl rt p\ 

/ {p + a)e^^^{x,t)dx + i{i + i) / / n^e^-^\e^\^dxdT + si / / {x + e)ps{e)e^dxdT 

Jo Jo Jo Jo Jo 

<Ci^T + Csl [ [ pOdxdT + Csl [ \\0\\]!L^. f p9^+^dxdT 
Jo Jo Jo Jo 

< {Q,T + CtI) + Csl ^ j^{p + a)9'+'dxdT. 
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Moreover, by using Gronwall's inequality, 



[\p + a)e^+\x, t)dx < {Ci^T + CtI) + {Ci^T + CTl)e^^^ 
Jo 



and 

On the other hand, by taking / — > cxd, we have 



\l^{Qt) ^ C't (3.28) 



where we have noted the fact 



1 

! + l 

:,T 

Moreover, by taking Z = 1 in the equation (|3.24|) . we obtain 



r {p + a)e\x,t)dx + \ f [\{Ki + K2\ps\^)\e^f + se^ps{e))dxdT <Ct, 

Jo ^ Jo Jo 

which imphes that 

\\Gx\\l2{Qj.), \\Pedx\\L2{Qj.) < Ct ■ (3.29) 

Secondly we present some estimates for p. By multiplying p on both sides of the equation 
()3.18p and integrating the resulting equation over Qt, with Gronwall's inequality we get 



sup 

o<t<r 



[ p'dx +11 \pA''dxdt < Cs,T + Cie, \\ipee,)e\\L^iQ^)) < Ce,T, (3.30) 

Jo Jo Jo 



which together with the Sobolev embedding inequality gives 

/ p'^dxdt < Ce,T- 
Jo 

Once again, multiplying p^ on both sides of the equation ()3.18p and integrating the resulting 
equation over Qt lead to 



sup / p'^dx + / / p'^\px\'^dxdt < Cs^T- (3.31) 
Ki<TJO Jo Jo 



™1 fl 

p'^dx+ / / 

0<t<TJo Jo Jo 

From ()3.28p . ()3.29p and ()3.30p . we conclude that {p, 9) is uniformly bounded in X'^ . Thus, by 
the Leray-Schauder fixed point theorem, there exists a fixed point {p^''^ , 6^''^) for the mapping 
M(-, •, 1) : X2 ^ X2 and {p'''' ,9"^'") is a solution of the system ([3:21) - ([33]) . 

3.3 Positivity of the approximate solutions 

Finally we prove the positivity of the approximate solutions {p^''^, 9'^'"). Let 9^ = 9e^ — 5. Then 
9^ is the solution of the following problem, 

{p + a)9f - {k'9% - [{e + {p9),)px + p{pe9.)e\ 9^ -{p + a)¥ - px%V9)9' 

+ q{9e\ 5)¥ = px%^)9e' + Xpx%V9)e' + (p + a)6 - (A + e-'5)ps{e-'6)e\ (3.32) 
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with the initial and boundary conditions 

' K^ei + = f3\9^e^ - 6), at x = 1, 

-K'ei + P"9^ = P\e^e^ -6), at x = 0, (3.33) 

_ e\x,0) = Ooeix) - S, for xen, 



where 



^ (A + e-'~e)ps{e-'~e) - (A + e-'5)ps{e-H) ^ ^ 
6-5 



By the assumption (|1.9|) . the right hand side of the equations (|3.32|) - (j3.33|) are nonnegative if 
5 is small enough (independent of e and v). Multiplying {0^)~ /{p + cr) on both sides of the 
equation (j3.32p and integrating the resulting equation over Qt, we derive ¥ > 0, i.e. B > e-'^6, 
which together with p.28p implies that 

Orp < e{x,t) < 9t for (x,t) G Qt- (3.34) 

where Orp and 9t are positive constants independent of e and v. 

For p, we define = p — 5. Then is the solution of the following equation 

pI - {{e + {pe),)pi), - {p\pe0.)e). + p'x'iVe) = x%Psi9)) + 6[{pee,),i - 5x'{Ve), (3.35) 

with the initial and boundary conditions 

{e + {p6)^)pi + p^{p,9^)e + a^p^ = a\p^ -5)-6{p,e^)„ at x = 1, 
-{e + {pe),)pi-p^{pee^)e + a^p^ = a\ffi -6) + 6{pee^,)e, at x = 0, (3.36) 
p^{x, 0) = poeix) — 6, for x S 0. 

Since x^(Ps(^)) > the right hand side of the equations (I3.35p - (I3.36[) are nonnegative if 

= mm< — , 



2' l + 2\\{ps9x)£\\cicQ^) 
in which case p^ > 0, or eqivalently p > 6. On the other hand, from (|3.29p we have 

Thus, there exists a positive constant p^ ^ such that 

P>P^^ for {x,t)eQT. (3.37) 



4 Global existence 



We have constructed an approximate solution {p^''^, 9"^'^) to the system (j3.ip and (|3.3p (or equvi- 
lently (|3.2p - (j3.3p ) in the last section. In this section, we prove the global existence of weak 
solutions for the system (ll.5p - ()1.8p . Firstly we fix e > and study the convergence as i/ — > 0. 

Since the system ()3.18p - ()3.19p reduces to (|3.2p - ()3.3p when s = 1, the uniform estimates 
()3.28p . (I3.29p . ()3.30p and ()3.34p given in the last section still hold for the approximate solution 
{p'^''^ ,9'^''^). We rewrite the first equation in (j3.2p by 

Pt = - fx + 9 
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with g uniformly bounded in L'^{Qt) and 

f = {e+ {p9)u)px + p{PeOx)e- 

Since p is uniformly bounded in L°° (I ; L'^ (0,)) n L'^{I\H^{Q)) ^ L^{Qt), we derive that 

\\piPe0x)s\\mQT)^ \\ip0)u\\Le(QT)^ \\ip0)upx\\ ^5 ^^^^ < Ce,T 

and 



From the first equation in (j3.ip we derive that 

II + 11^5/4^^.^^-1,5/4 (^^^ < C,,T 

where we have noted p.28p . and moreover, from (j3.30p . we observe that p^''^ is uniformly bounded 
in L^{I-L^{Q)) n L'^{r,H\n)) and pf" is uniformly bounded in L5/4(j. p^-i.5/4^^^^_ -^^.^^ 
Aubin-Lions lemma, we conclude that there exists a sequence Vj — > such that 

pe,u, _^ pS strongly in LP{Qt) (VI < p < 6), 
p^,'^. ^ strongly in L'^ {I; 0(11)), 

(4-1) 

p^^'^j^p^ weakly in L^{I,H^{n)), 
p't'^'^pt weakly in L5/4(/; VFo"^'^/^(17)) 

and 

p^:-.(0,.) ^p^(0,.) and p^-'^^Xl,-) strongly in L^{0,T). 

Similarly, by noting the uniform estimates (|3.28p . (j3.29p and (j3.34p . we conclude that there 
exists a subsequence of 0^'^^ (also denoted by 9^''^^) such that 

Qe,uj _^ pS strongly in WIQt) (VI < p < oo), 
Qe,uj^Qe strongly in L2(/;C(n)), 

(4.2) 

0e,yj^e^ weakly in L^{I,H^{n)), 

(p^.^-j^''^^ +(7r''^^)t^(p"r + weakly in L^/^{I;W~^'^^\n)) 

and 

r'^^ (0,-) ^ 6*^(0, •) and ^'^^ (1, •) ^ 6*^(1, •) strongly in 1^(0, T), 1 < p < oo. 
Since {p^''^^ ,0^''^^) is a strong solution of the system (13. ip and p.3p . it satisfies 

a^ip''''^{0,t) - p'^)(P{0,t)dt + [ a\p''''^{l,t) - p^)^il,t)dt 



+ r f pY'ct>dxdt+ r [{e + {p''-^e'^^%)p:/' +p''''^{pr'e7')e](kxdxdt 
Jo Jn Jo 

[ x'iPs{9'''''))dxdt- [ [ p^^''^x%Ve^)(l)dxdt 
Jn Jo Jn 
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and 
rT rl 



^0 



T 



+ / a"{p''^''^{0,t)- p")e'^''^{0,t)i^{0,t)dt+ / a\p'^''^il,t) - p')9'^''^{l,t)ij{l,t)dt 
Jo Jo 

+ 11 K'e7'^.dxdt+ [ [ [{e + {p''''^e''''%)p'/^e''''^ +p''''^{p7^e7')e0''''']^xdxdt 



+ 



JO Jo Jo 

[ [ {X + e^^''^)psiO''''^)'^]dxdt 
Jo Jo 



(4.4) 



= x[ [ p''^''^x'{V¥^)ijdxdt + X [ [ e^^''^x%PsiO''''')Wxdt, 
Jo Jo Jo Jo 

for any L^{I;W^'^{Q)). By taking the limit j oo, we obtain a global weak solution 

{p'^,6^) to the approximate system 

Pt - ((e + p6)px)x - {p{PeOx)e)x = -^e, 

{pO + ae)t - {k^Ox), - ((e + p9))pxe)x - {p{pe9cc)ee)x (4.3) 

= XT, + iX + 9) ix%Psm-Ps{e)), 

with the boundary and initial conditions 

{e + p9)px + P(p^x)4=i = "'(Z'' - /'(I' *))' 
{e + p9)px + p{p9x)e\,^Q = a°(p(0, t) - f), 
p{x,0) = poeix) := po * r]e(x) +e, 
k'9x\x=i = P\9' -9{l,t)), 
k'9x\x=o = P''{9{0,t)-9'>), 
9ix,0) = 9o,ix) ■.= 9o*r],{x). 

Secondly, we study the convergence as e ^ 0. To take the limit e — > 0, we need more 
uniform estimates for p with respect to e. 

Clearly the system ()3.2p - ()3.3p reduces to the system ()4.3p - (|4.4p when u = 0. Then the 
uniform estimates p.l6p and ()3.17p hold for the obtained solution {p^,9^). From ()3.3ip we see 
that 

WpGPxWl^Qt) ^ Ce,T 

and from the first equation of ()4.3p we deduce that pt £ L^{I; HQ^{il.)). Note that Inp G 
L^{I; H^{Q,)). By multiplying the first equation of (j4.3p by Inp and integrating the equation 
over Qt, we arrive at 

/ plnp{x,t)dx— I p{x,t)dx+ f [epx + p9px + p{ps9x)e]^T'^ P dr + f f 9p^dxdT 
Jo Jo Jo Jo Jo 

< / pQelnpQeix)dx - / poedx- / / {pe9x)ePxdxdT - / / {pV9 - Ps{9)) In p dxdr. 
Jo Jo Jo Jo Jo Jo 
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Since t i t i t i 2 

\{Pee^)epx\dxdT <\f f 0pldxdT + l [ [ K^^^l^rfxdr 



'0 -'0 Jo Jo ^ Jo Jo 



< 

- 2 



1 



2 

we get 



< ^ / / Opldxdr + Ct ■ 
Jo 



f pin p{x,t)dx -\ — / / 9p'^dxdT + / / plnpdxdr 

J[o,i]n{p>i} Jo Jo J ^[o,i]x[o,i]n{p>i} 

< / poe\ln pos\ix)dx + 
Jo Jk 



p\ In p\ {x, t)dx 

[o,i]n{p<i} 



+ / / p\liip\dxdT+ / / ps{6)liipdxdT + Ct 

J i[o,i]x[o,t]n{p<i} J ^[o,i]x[o,t]n{p>i} 

<Ct, 

which, together with ()3.34p . leads to 

||plnp||ioo(o,T;Li(n)), \\Px\\l2{Qt) < Ct ■ (4.5) 
From the inequahties (13.27P and (14. 5p we derive 

\\p\\L^o,T;m{n)) < Ct (4.6) 

and 

iip|iio=(n)<iipiiii(f.) + iiHiK(n)ii/'-iiL4) 

1, 



which results in 



Moreover, we have 



<CT + ^\\p\\U^n) + CT\\Px\\h^n)^ 



£ \\p\\looin)dt <Ct + Ct £ \\p.\\h^n)dt < Ct- 



I I p'^dxdt<(^l Qup^^ pdx^<CT, 



(4.7) 



i.e. p is uniformly bounded in L'^(Qt)- 
Finally, we let 

Then Bi B2 ^ B^ and {p^} is uniformly bounded in L^(/; B2) H L^{I; Bi). From the first 

equation in ()3.2p . i.e. 



Pt 



epx + p0Px + Pip0x)e] - Px'iVe) + x%Ps{e)) , 

- X 
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we observe that {pf} is uniformly bounded in L^/^{I; B3). By Aubin-Lions lemma, {/o^} is 
relatively compact in for {1 < p < 4). Thus, there exists a sequence p^^ such that 

limj^oo £j = and 

p^3 p strongly in L'^{Q)) (VI < p < 4), 

p^^ ^ p strongly in L^(I, (7(17)), 

(4.8) 

p^^^p weakly in L^{I,H^{n)), 
pl'^Pt weakly in L^/^^{I-W~^'^^^{Vt)). 
Similarly, by (|3.16p and (|3.17p . there exists a subsequence of 9^^ (also denoted by 6^^) such that 

9 strongly in LP{Qt) (VI < p < 00), 
r^^e strongly in L^{I,C{Q)), 

(4.9) 

9^3 ^9 weakly in L^{I,H^{Q)), ^ ' 

(p^^r^ +(jr^)t^(p0 + CT0)t weakly in L^I^{1;W~^''^I'\Vl)). 

Now we take the limit j — > 00 and by (14. Sp and (j4.9l) . we obtain a weak solution (p, ^) which 
satisfies ([23]) and (12:2]) . 

Acknowledgements The authors wish to thank Professors P. Lei, T. Yang, G. Yuan and 
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